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Abstract: Recently we developed a new framework in Hirz et al. (2015 1 to model stochastic mortality us¬ 
ing extended CreditRisk+ methodology which is very different from traditional time series methods used for 
mortality modelling previously. In this framework, deaths are driven by common latent stochastic risk factors 
which may be interpreted as death causes like neoplasms, circulatory diseases or idiosyncratic components. 
These common factors introduce dependence between policyholders in the annuity portfolios or between death 
events in population. This framework can be used to construct life tables based on mortality rate forecast. It 
also provides an efficient, numerically stable algorithm for an exact calculation of the one-period loss distri¬ 
bution of annuities or life insurance products portfolios and associated risk measures such as value-at-risk and 
expected shortfall required by many regulators. Moreover this framework allows stress testing and, therefore, 
offers insight into how certain health scenarios influence annuity payments of an insurer. Such scenarios may 
include improvement in health treatments or better medication. In this paper, using publicly available data for 
Australia, we estimate the model using Markov chain Monte Carlo method to identify leading death causes 
across all age groups including long term forecast for 2031 and 2051. On top of general reduced mortality, 
the proportion of deaths for certain certain causes has changed massively over the period 1987 to 2011. Our 
model forecasts suggest that if these trends persist, then the future gives a whole new picture of mortality for 
people aged above 40 years. Neoplasms will become the overall number-one death cause. Moreover, deaths 
due to mental and behavioural disorders are very likely to surge whilst deaths due to circulatory diseases will 
tend to decrease. This potential increase in deaths due to mental and behavioural disorders for older ages will 
have a massive impact on social systems as, typically, such patients need long-term geriatric care. 


Keywords: Extended CreditRisk+, stochastic mortality model, life tables, annuity portfolios, life insurance 
portfolios, longevity risk, risk management, estimation of extended CreditRisk+, Markov chain 
Monte Carlo. 
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1 INTRODUCTION 


Mortality modelling has a very long history. Numerous deterministic survival models for mortality intensity 
have been developed since Benjamin Gompertz suggested an exponential increase in death rates with age 
iGompertz’ law of mortality) in 1825. However, stochastic modelling of mortality is more modern develop¬ 


ment over last twenty years. One of the standard benchmarks a model developed by Lee and Carter (1992 1 , for 
modelling log death rates of people across different ages affected by one time dependent latent factor. There 


are many extensions introducing several factors and special treatment for cohort effects; for review, see Booth 
[and Tickl^ ( |2008| l. This stochastic modelling has become increasingly important in the financial industry and 
government departments since it is observed over the last several decades that life expectancy is typically un¬ 
derestimated. Moreover, new regulatory requirements such as Basel III and Solvency II and the financial crises 
of 2008 also call for stochastic mortality modelling. 

Life insurers and pension funds usually use deterministic first-order life tables to derive premiums, forecasts, 
risk measures for portfolios and other related quantities. These first-order life tables are derived from second- 
order life tables (best estimates of the current mortality of a population) plus artificially added risk margins 
associated with longevity, size of the company, selection phenomena, estimation and various other sources, 
see, for example, Pasdika and Wolff ( 2005| l. The risk margins described there often lack stochastic foundation 
and are certainly not consistently appropriate for all companies due to a possibly twisted mix of these risks. 
Moreover, drastic shifts in death rates due to certain death causes are observed over the past decades. This 
phenomenon is usually not captured by generation life tables which incorporate only an overall trend in death 
probabilities. As an illustration of this fact, Figure[2shows death rates based on Australian data. 
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Figure I. Australian death rates for mental and behavioural disorders (left), as well as for circulatory diseases 
(right) from 1987 to 2011 for age categories 70-79 years, and 60-69 years with both genders. 


In our recent paper Hirz et al. (2015 i, a collective risk model called extended CreditRisk+ was adapted to 
modelling mortality and calculation of loss distribution for annuity portfolios. This model setup is very differ¬ 
ent from traditional time series models developed for mortality modelling previously. As the name suggests, 
it is a credit risk model used to derive loss distributions of credit portfolios and originates from the classical 
CreditRisk+ model which was introduced by Credit Suisse First Boston (1997 1 . Within credit risk models it 
is classified as a Poisson mixture model. Identifying default with death makes the model perfectly applicable 
for various kinds of life insurance portfolios and annuity portfolios. 


In this paper, using publicly available data for Australia, we estimate model parameters via a Markov chain 
Monte Carlo (MCMC) method to identify leading death causes across all age groups including long term fore¬ 
cast for 2031 and 2051. Sectionj^presents a model formulation. Likelihood function and MCMC method for 
estimation are summarised in Sectionj^ Results are presented in Section|^and Section|^contains concluding 
remarks. 
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2 Model 


We follow the model developed in |Hirz et al.| ( |2015] l where the reader can find precise mathematical formula¬ 
tion. Here we provide a brief outline of the main model blocks and assumptions. 


Let m} denote the set of people in population (or policyholders in the annuity portfolio) with death 

indicators Ni,, Nm over some period of time (i.e. one year if we model annual death rates or annual 
portfolio loss distribution). Event Ni — 0 indicates no death for i G m}. In reality, death indicators 

Ni are Bernoulli variables as each person can die only once. However for easier calibration and subsequent 
calculations a Poisson approximation is assumed that does not lead to material difference in the case of large 
number of people in a group as death probabilities over one year are small, see discussions and numerical 
examples in Hirz et al. (2015|l. 


Denote latent stochastic risk factors identified with causes of death (such as neoplasms, cardiovascular diseases 
or idiosyncratic components) by A^, fc = 1,... ,K which are modelled as independent random variables from 
gamma distribution with mean 1 and variances a\,k = 1,..., AT. These stochastic risk factors are designed to 
model effects which simultaneously influence death probabilities of many people due to a common exposure 
to the same type of risk. Risk index fc = 0 is reserved to represent idiosyncratic risk. 


Denote corresponding weights (vulnerability of policyholder i to risk factor k) as Wi fe and by model construc¬ 
tion Wifi -f • • • -f Wi^K = 1- In average, when we consider homogeneous groups of people this weight is the 
fraction of people dying from death cause k compared to all deaths in the group. 

Eor every person i G to}, the total number of deaths Ni is split up additively according to risk factors 

as 


Ni — Nifi -f • • • -f Ni^ii, (1) 

i.e. A^i fe is the number of deaths of a person i due to risk factor k. It is assumed that death indicators for 
idiosyncratic risk iVi o,..., N^.o are independent from one another, as well as all other random variables and 
Nifi is Poisson distributed with intensity qiWifi, where qi denotes the probability of death of a person i. 

Conditionally (given risk factors), death indicators Ni^k are independent and Poisson distributed with random 
intensity qiWi^k^k, i-S-, 


/ m K \ m K 

(n Ai,...,Aif j=]^]^ 


^ — qiWi^k 


Afc {qrWi.k^kY 


nik'- 


( 2 ) 


Thus, by our model construction, E[A^i] = qi{wifi -f • • • -f Wi^K) = qi- Note that under this model, loss 
distribution for portfolio of annuities can be calculated exactly and efficiently using Panjer recursion algorithm 
instead of approximate and slow Monte Carlo method; for more details, see Hirz et al. (2015|l. 


Consider time periods t G {1,... ,T,...} where T corresponds to the last year in the dataset. To model trends 
in time, death probabilities qi and weights Wi fe are assumed to be time dependent. If the modeller is interested 
in prediction over one year, then simple linear functions of time can be appropriate, however for long term 
forecasting these are designed more carefully as 


q^t) = [a, + YiTt^.^riiit)), w,^k{t) = ' 


exp{u^^k + Vi,kTpk,1|>kit)) 


(3) 


with model parameters ai, fii: G G ® ^nd rji G (0, oo), as well as Uifi^ Vifi, (j>o, ..., Ui^x, Vi^x, 4>k C R, and 
■00,..., € (Oj oo) to be estimated using data. Here, is Laplace distribution and is trend 

reduction function formally defined as 

i-f ^sign(a:)(l - exp(-|x|)) , x gM., (4) 




- arctan(^ + rjt), 


f e K, (C, ? 7 ) C K X (0, oo). 


(5) 


Eor a; < 0, Laplace distribution (|4| becomes exp(a;)/2 and thus corresponds to modelling of log death rates. 
Trend reduction function is motivated by Kainhofer et al. ( 2006[ section 4.6.2) where linear time trend was 
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replaced by time shift 7o,,;(f) with 77 = Then, parameter rj gives the inverse of the time to when an 
initial trend is halved. Parameter ( on the other hand gives the shift on the arctangent curve. Note that 
lima;_>±oo arctan(a:) = ±f and thus death probabilities and weights ([^ are non-degenerate, i.e., do not hit 
zero or one, as f — 00 . 

3 Data Likelihood and MCMC estimation 


For every age category a G {1,..., A}, gender g G {f, m} and year t G {1,..., T} with T > 2 the database 
is assumed to contain historical population counts ma,g{t) and historical number of deaths na,g,k{t) due to 
underlying death cause k G {0, 1,... ,K}. An underlying death cause is to be understood as the disease or 
injury that initiated the train of morbid events leading directly to death. 

The observations of historical annual deaths with age o G A}, gender g G {f, m}, due to 

death cause k G {0, ... ,K} and at time f G {1,..., T} correspond to realisations of the random variable 

Na,g,k{t)-- ^ 

where Ma^g{t) C m{t)} denotes the set of people (policyholders) of specified age group and gender. 

Note that ^(f) is the number of deaths of a person i due to death cause k in year t. Death cause zero 
corresponds to ill-defined and not elsewhere reported deaths, i.e., idiosyncratic components, but could also be 
put as another stochastic risk factor. 

It is assumed that for alH G {1,..., T}, quantities qi{t) and corresponding weights Wi^k{t), respectively, are 
the same for all people i G rn{t)} within the same age category a, same gender g and with respect to 

the same risk factor Ak{t). Therefore, we can set qa,g{t) := qi{t) and Wa,g,k{t) ■= Wi^kit) for a representative 
policyholder i of age category a and gender g with respect to risk factor Afc(f). Define 

.4 

a—l {f ,ni} 

as well as Pa,g,k{t) := ma,g{t)qa,g{t)wa,g,k{t) and 

.4 

Pk{t):='y2 

a—l {f ,ni} 


Then, the likelihood function £{n\9g, 6 ^,( 7 ) of parameters Og := {a, /3, C, v), as well as 6 ^ ■= (u, v, (j>, ip) 
and (T := (cti, ..., aK) given mortality data n := {na,g,k{t), a — 1,... ,A,g G {f, m}, t = 1,..., T) can 
be easily calculated. Given risk factors Ai,..., Ak the conditional likelihood is just a product of appropriate 
Poisson probabilities, then the unconditional likelihood can be found by integrating out the gamma distributed 
independent Ai,..., Afc and is given by 


(.{n\eg,e^,cT) = wi(W 

t=l \ a—l pG{f,ni} 


K 

n 


'^a,g,k (f) ■ 

r{l/al+nk{t)) 




(t) 


A 


xn n 

a—l g^ {f ,ni} 


^0,3,/ J 


(6) 


It is also possible to derive closed form maximum a posteriori (MAP) estimates for risk factor realisations and 
their variances; for details of MAP estimates and likelihood derivation, see |Hirz et al.| ( |20f5] l. 

Once the likelihood is calculated in closed form, it is straightforward to apply MCMC method to get many 
samples 0* = (0*, cr®), i = 1,2,... from the posterior density 7 r( 0 |n) oc £{n\0)TT{0). The mean over 

these samples provides good point estimate for the model parameters and sampled posterior density can be 
used to estimate parameter uncertainty. We use non-informative (uniform) prior density Tr{0) for model pa¬ 
rameters; in this case the mode of the posterior samples corresponds to a maximum likelihood estimate and 
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inference is based on data with non-material impact from the prior. There are many MCMC algorithms that 
can be applied. In this study we implement a well known random walk Metropolis-Hastings within Gibbs 
algorithm as described e.g. in ( |Cruz et [2015 section 7.4.4). Also, we utilize truncated normal distributions 
as proposal distributions. The method requires a certain burn-in period until the generated chain becomes 
stationary. Estimates derived by matching of moments as described in |Hirz et al.| ( |2015| l can be used as initial 
values to ensure a shorter burn-in period. To reduce long computational times, one can run several independent 
MCMC chains with different starting points on different CPUs in a parallel way. We have tested our imple¬ 
mentation via simulated experiments where we simulate risk factors and data from the model with known 
parameters and then estimate risk factors and parameters. 


4 Study of Australian mortality data 


We applied the above described model and MCMC estimation procedure to Australian death rate data for the 
period 1987 to 2011. Data source for historical Australian population, categorised by age and gender, is taken 
from the website of Australian Bureau of Statistic http;//www.abs.gov.au and data for the number of deaths 
categorised by death cause and divided into nine age categories (0-9 years, 10-19 years, 20-29 years, 30-39 
years, 40—49 years, 50-59 years, 60-69 years, 70-79 years and 80 h- years, denoted by oi,..., og, respectively) 
for each gender is taken from the Australian Institute of Health and Welfare website http://www.aihw.gov.au 


Data handling needs some care as there was a change in International Classification of Diseases in Australia 
in 1997 as explained at the website of the ABS. As a result, for the period 1987 to 1996, death counts have 
to be multiplied by corresponding comparability factors and rounded to the nearest integer in order to avoid 
data inconsistencies. The provided death data is divided into 19 different death causes where we identify the 
following ten of them with common non-idiosyncratic risk factors (comparability factor is given in brackets): 
certain infectious and parasitic diseases (1.25), neoplasms (1.0), endocrine, nutritional and metabolic diseases 
(1.01), mental and behavioural disorders (0.78), diseases of the nervous system (1.2), circulatory diseases 
(1.0), diseases of the respiratory system (0.91), diseases of the digestive system (1.05), external causes of 
injury and poisoning (1.06), diseases of the genitourinary system (1.14). We merge the remaining eight death 
causes to idiosyncratic risk, termed as not elsewhere defined, as their individual contributions to overall death 
counts are small for all categories. 


Trend reduction parameters are fixed a priori with values C,ai,g = fik = 0, as well as rja-^g = V'fe = i|o 
ages i G 9}, g G {f, m} and k G {0,..., K} with K = 10. Thus, we have to estimate 442 parameters, 

36 of which can be chosen arbitrarily as the system is overdetermined. Eixing the parameter for trend reduction 
makes estimation more stable and does not influence results significantly for mid-term forecasts. The value 
of for rjai,g and ipk is an average approximation to trend reduction observed in Australia which is usually 
higher for very old and mid ages, see Hirz et al. (20151. Based on 35 000 MCMC steps with burn-in period 
of 5 000 we are able to derive estimates of all parameters. The results of estimation for leading death cause 
weightings including forecast for 2031 and 2051 across all age groups, using ([^ and MCMC samples of the 
parameters in the usual way, are presented in Tables [T] and for males and females, respectively. On top of 
general reduced mortality, the proportion of deaths for certain certain causes has changed massively over the 
period 1987 to 2011. Our model forecasts suggest that if these trends persist, then the future gives a whole 
new picture of mortality for older ages. Weightings for death causes up to the age of 39 are estimated to be 
relatively stable over time. Note that for the youngest age group, the leading death cause is idiosyncratic risk 
since it contains deaths related to the perinatal period, as well as all congenital malformations, deformations 
and chromosomal abnormalities. Eor ages above 40 years we derive an overall relative decrease in deaths due 
to circulatory diseases with small estimation error (see MCMC quantiles in Tables [^and]^. Deaths due to 
neoplasms will become the overall number-one death cause, but weights are relatively stable. Weightings for 
deaths due to mental and behavioural disorders show the most dramatic change and this cause will most likely 
become one of the leading death cause for very old people. This potential increase in deaths due to mental and 
behavioural disorders for older ages will have a massive impact on social systems as, typically, such patients 
need long-term geriatric care. But it should be outlined that estimation errors are high for this risk factor as 
can be seen in the wide confidence intervals in Tables [T] and |2] 


As outlined in Hirz et al. (20151 using similar data, model validation techniques strongly suggest that our 
proposed model is suitable for describing the observed data. These validation techniques include a test for 
checking whether sample covariances amongst deaths of the same cause lie within confidence intervals, a 
test for independence amongst death counts of different death causes using a f-test, as well as tests for serial 
correlation and for distributional assumptions on risk factors. 
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5 CONCLUSIONS 

This paper presents estimation results of leading death causes for all age groups in Australia including fore¬ 
cast for 2031 and 2051. The adopted modelling framework of extended CreditRisk+ is very flexible. In a 
general form, extended CreditRisk+ risk factors can be used to model risk groups with simultaneous deaths of 
policyholders in the group (e.g. a couple dying in a car crash, people living near a volcano, virus outbreaks), 
moreover dependence (negative and positive) between death causes can be introduced. Development of this 
general case for mortality modelling is a subject of further research. 
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Table 1. Leading death causes with weights estimated by MCMC sample average with 5 and 95 percent 
confidence intervals of MCMC posterior in brackets for males of all age categories in 2011, 2031 and 2051. 



2011 

2031 

2051 


1. 

not elsewhere: 0.669 (o g5g) 

not elsewhere: 0.593 (g'ffg) 

not elsewhere: 0.456 (g'lfg) 

0-9 years 

2. 

external: 0.108 (g'JJf) 

circulatory: 0.077 (g'Jfg) 

circulatory: 0.157 


3. 

neoplasms: 0.053 (o qI®) 

external: 0.074 (0;0®|) 

digestive: 0.133 (o:f||) 


1. 

external: 0.651 (o gf®) 

external: 0.544 (f m) 

external: 0.426 (g ||g) 

10-19 years 

2. 

neoplasms: 0.099 (g jgo) 

nervous: 0.117 (g'J|f) 

nervous: 0.155 (g fjl) 


3. 

nervous: 0.079 (g off) 

neoplasms: 0.110 (g'J|g) 

neoplasms: 0.113 (g gyg) 


1. 

external: 0.732 (o'72i) 

external: 0.701 (0 672) 

external: 0.653 (g fgg) 

20-29 years 

2. 

neoplasms: 0.072 (q qIy) 

neoplasms: 0.084 (0 072) 

neoplasms: 0.094 (o J||) 


3. 

circulatory: 0.054 (0 Q49) 

circulatory: 0.072 (0 059) 

circulatory: 0.093 (g jfg) 


1. 

external: 0.590 (g'gyg) 

external: 0.648 (g'fgf) 

external: 0.679 (g glf) 

30-39 years 

2. 

circulatory: 0.125 (g jgo) 

circulatory: 0.125 (g'}|f) 

circulatory: 0.120 (o ||g) 


3. 

neoplasms: 0.119 (g jfl) 

neoplasms: 0.095 (g'Jgf) 

neoplasms: 0.074 (g ggg) 


1. 

external: 0.307 (g'lgl) 

external: 0.359 (g'337) 

external: 0.398 (g'lff) 

40^9 years 

2. 

neoplasms: 0.256 (g'^tg) 

neoplasms: 0.233 (g'||g) 

neoplasms: 0.204 (g f||) 


3. 

circulatory: 0.205 (g igg) 

circulatory: 0.148 (g'Hf) 

circulatory: 0.104 (g jgf) 


1. 

neoplasms: 0.428 (o'lfs) 

neoplasms: 0.437 (g'|ff) 

neoplasms: 0.409 (g's^g) 

50-59 years 

2. 

circulatory: 0.229 

external: 0.138 (g'||g) 

external: 0.148 (q'i 3 q) 


3. 

external: 0.117 (g Jis) 

circulatory: 0.135 (g'}||) 

infectious: 0.102 (q'Jq®) 


1. 

neoplasms: 0.500 (g lgg) 

neoplasms: 0.568 (g'ffg) 

neoplasms: 0.583 (g fgl) 

60-69 years 

2. 

circulatory: 0.233 (g'^gg) 

circulatory: 0.118 (g'}||) 

endocrine: 0.088 (g jff) 


3. 

respiratory: 0.061 (g gfg) 

endocrine: 0.067 (g'ggf) 

nervous: 0.067 (g gfg) 


1. 

neoplasms: 0.432 (g igs) 

neoplasms: 0.514 (g'fgl) 

neoplasms: 0.538 (g ffl) 

70-79 years 

2. 

circulatory: 0.270 (g gg|) 

circulatory: 0.133 (g'Jfg) 

endocrine: 0.109 (g jgl) 


3. 

respiratory: 0.092 (g ggg) 

endocrine: 0.078 (g'gfl) 

nervous: 0.077 (g ggg) 


1. 

circulatory: 0.376 (g pl) 

neoplasms: 0.297 (g'ggf) 

neoplasms: 0.301 (g pf) 

80+ years 

2. 

neoplasms: 0.259 (g gfg) 

circulatory: 0.244 (g'lff) 

mental: 0.160 ([j;f22) 


3. 

respiratory: 0.112 (g }g|) 

mental: 0.091 (g gH) 

circulatory: 0.146 (g lff) 
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Table 2. Leading death causes with weights estimated by MCMC sample average with 5 and 95 percent 
confidence intervals of MCMC posterior in brackets for females of all age categories in 2011, 2031 and 2051. 



2011 

2031 

2051 


1. 

not elsewhere: 0.685 (g g®®) 

not elsewhere: 0.631 (g'lgf) 

not elsewhere: 0.531 (Q' 44 g) 

0-9 years 

2. 

external: 0.090 (g gll) 

circulatory: 0.086 (g'ggf) 

circulatory: 0.184 (g fjf) 


3. 

neoplasms: 0.057 (g gfl) 

neoplasms: 0.070 (g glf) 

neoplasms: 0.076 (® ggg) 


1. 

external: 0.507 (g lgg) 

external: 0.439 (g'Hg) 

external: 0.363 (g lgg) 

10-19 years 

2. 

neoplasms: 0.152 (g }®®) 

neoplasms: 0.169 (g'fgf) 

neoplasms: 0.175 (g ff®) 


3. 

not elsewhere: 0.084 (g g®®) 

nervous: 0.104 (g'gfg) 

nervous: 0.125 (® g®g) 


1. 

external: 0.508 (g igf) 

external: 0.454 (g'|f®) 

external: 0.397 (g HI) 

20-29 years 

2. 

neoplasms: 0.167 (g jgg) 

neoplasms: 0.185 (g'fg®) 

neoplasms: 0.196 (g fl®) 


3. 

circulatory: 0.081 (g g®®) 

circulatory: 0.102 (g ggg) 

circulatory: 0.122 (g ggg) 


1. 

external: 0.337 (g'H®) 

external: 0.366 (0 333 ) 

external: 0.381 (g'glg) 

30-39 years 

2. 

neoplasms: 0.297 (g |gg) 

neoplasms: 0.233 (0 212 ) 

neoplasms: 0.180 (g f®!) 


3. 

circulatory: 0.120 (g }f|) 

circulatory: 0.124 (q'}q 7 ) 

circulatory: 0.123 (g jp) 


1. 

neoplasms: 0.480 (g ifl) 

neoplasms: 0.411 (g'|gg) 

neoplasms: 0.338 (g'lgg) 

40^9 years 

2. 

external: 0.164 (g jg^j.) 

external: 0.189 (g'fgg) 

external: 0.204 


3. 

circulatory: 0.133 (g }®?) 

circulatory: 0.113 (g'}g|) 

circulatory: 0.093 (g j®!) 


1. 

neoplasms: 0.588 (g fgx) 

neoplasms: 0.586 (g'®gg) 

neoplasms: 0.553 (g'lgg) 

50-59 years 

2. 

circulatory: 0.120 (g jxg) 

external: 0.092 (g'Jgg) 

external: 0.105 (® J||) 


3. 

external: 0.076 (g g^l) 

circulatory: 0.063 (g g®!) 

not elsewhere: 0.063 (g glg) 


1. 

neoplasms: 0.566 (g fgi) 

neoplasms: 0.640 (g ggg) 

neoplasms: 0.666 (® g|g) 

60-69 years 

2. 

circulatory: 0.157 (g jfg) 

respiratory: 0.075 (g;ggg) 

nervous: 0.067 (g ®®!) 


3. 

respiratory: 0.079 (g gfg) 

circulatory: 0.066 (g'ggl) 

respiratory: 0.067 (g gg^) 


1. 

neoplasms: 0.405 (g |gg) 

neoplasms: 0.478 (g;!!®) 

neoplasms: 0.495 (g lgg) 

70-79 years 

2. 

circulatory: 0.256 (g' 252 ) 

respiratory: 0.114 (g'Jgg) 

respiratory: 0.115 (® (q®) 


3. 

respiratory: 0.100 (g'Jg|) 

circulatory: 0.111 (g'JJg) 

nervous: 0.094 (g jgf) 


1. 

circulatory: 0.430 (g l®!) 

circulatory: 0.252 (o'^sg) 

mental: 0.255 (® fgg) 

80+ years 

2. 

neoplasms: 0.167 (g }®|) 

neoplasms: 0.180 (o' 170 ) 

neoplasms: 0.163 (0 ( 43 ) 


3. 

respiratory: 0.093 (g ggg) 

mental: 0.151 (o'Jll) 

circulatory: 0.130 (® (f|) 
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